Resonances in a circular dielectric cavity 
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Abstract 

We study resonance distributions in a circular dielectric cavity. It is shown that the decay-rate 
distribution has a peak structure and the details of the peak are consistent with the classical 
survival probability time distribution. We also investigate the behavior of the complex resonance 
positions at the small opening limit (n — > oo, n is the refractive index of the cavity). At the large 
n limit, the real part of complex resonance positions approaches the solutions with different m of 
Dirichlet problem with a scale n~ 2 and the imaginary part goes zero as n~ 2m for TM and n~ 2 ( m+1 ) 
for TE polarization, where m is the order of the resonance. 
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Two-dimensional (2-D) billiards with various boundary geometries have been widely stud- 
ied in quantum chaos community because of the ease in analysis and rich interesting phe- 
nomena such as classical and wave chaos and wave localization etc. There are many relevant 
experimental realizations, e.g., surface waves microwave billiards , mesoscopic struc- 



tures 



0,flB 

, and semiconductor microcavities , [h], • In experiments, the systems 



measured are inevitably coupled with environment, i.e., they are open. Even in the case of 
a small coupling, it is sometimes difficult to explain the open system based on the physical 
properties of the corresponding closed system without a clear understanding of the corre- 
spondence between them. 

The decay-rate statistics is determined by interpaly of openness and dynamics of the 
system concerned 12J. In practical purpose, it is, however, not a simple task to obtain many 
decay rates enough to apply statistical processes for a chaotic open quantum system, such as 
a chaotic dielectric cavity, because of heavy numerical tasks. From the viewpoint of a easy 
analysis, the circular dielectric cavity is a good open system, due to its simple geometry, 
to study the statistics of decay rates which originate from the dielectric property and the 
relationship to the corresponding closed billiard with the Dirichlet boundary condition. In 
dielectric cavities the degree of opening is determined by the refractive index n and the small 
opening limit corresponds to n — > oo. 

In this letter, we answer the following two questions about resonances in a circular di- 
electric cavity: (i) How are the imaginary values of resonance positions, which represent 
the decay rates of the resonances, in the circular dielectric cavity distributed? (ii) Would 
its resonance positions approach the eigenvalues of the corresponding billiard at the small 
opening limit? For the first question, we obtain the complex resonance positions and show 
that the distribution of imaginary values is consistent with survival probability time distri- 
bution (SPTD) which represents the decay property due to classical ray escapes [131 ] . As for 
the second question, we show that the resonance of TM polarization with angular quantum 
number m in circular dielectric cavity approaches the eigenvalues with m — 1, not m, of 
Dirichlet problem at the small opening limit (n — > oo). 

From the mathematical viewpoint, the only difference between closed billiard and open 



dielectric cavity problems is the boundary condition applied to the Helmholtz equation 14| , 



(V 2 + ^ = 0, (1) 



2 



where k is the wavenumber inside the billiard or the dielectric cavity. In the closed billiard 
problem, the typical boundary conditions are Dirichlet and Neumann boundary conditions, 

(Dirichlet) ip(r b ) = 0, (2) 
(Neumann) d n ip(rb) = 0, 

where r& denotes the boundary vector and d n is the normal derivative to the boundary 
For an arbitrary boundary shape we can find real eigenvalues k n and the eigenfunctions 
ip n ( r ) satisfying the corresponding boundary condition. In the dielectric cavity case, we 
have to apply different boundary conditions depending on the polarization. For TM (TE) 
polarization, the electric (magnetic) field normal to the 2-D cavity and its normal derivative 
(its normal derivative divided n 2 ) are continuous at the boundary interface, i.e., 

out 

(TE) n *$™( T *) - ^out^outin) 



d n ipin(n) d n i{j ou t(r b ) 
where ^ n (r) and "i/w(r) are wave functions inside and outside the dielectric cavity, respec- 
tively, and the refractive indices of the cavity and environment are n in and n out . We will 
set rti n = n and n out = 1 throughout the letter. The resonance positions k, obtained from 
the above boundary conditions, are complex numbers with negative imaginary values. The 
quality (Q) factor of the corresponding resonance mode is defined as — Re(/c)/2Im(/c). High 
Q resonance modes, therefore, have smaller absolute values of Im(fc) and mean the resonance 
modes well confined inside the cavity. 

For the circular boundary geometry with a radius R, the Helmholtz equation, Eq.([T]), 
can separate into the angular and radial equations due to the rotational symmetry, and 
the solution of the angular equation is simply given by e tm< ^ (m is an integer). The radial 
equation is represented by Bessel differential equation [151 ]. and relevant solution inside both 
the circular billiard and dielectric cavity is Bessel function, J m (kr). In the circular dielectric 
cavity, the solution outside the cavity should be outgoing wave which is represented by 
Hankel function of the first kind, H^(kr). Using these solutions, the boundary conditions 



(Eq.(T5]), Eq.(j3])) become 

u no, m 



(Dirichlet) J m (k) = 0, (4) 
(Neumann) kJ m _i{k) — mJ m (k) = 0, 



(TM) nJ m _ x - J m {k)H^\{k ) = 0, 

(TE) nJ^H^lAh) - J m ^(k)H^(k ) 

= ^(n--)J m (k)HW(k ), 
Kq n 

where we take R = 1 without loss of generality, and k Q is the vacuum wavenumber and 
ko = k/n. In general, the separability of the Helmholtz equation implies that the circular 
system is integrable and there are good quantum numbers specifying the solution. In our 
case, the solution can be classified by the angular quantum number m and the radial quantum 
number I as h( m ,t)- Figure 1 shows eigenf unctions satisfying Dirichlet (a) and Neumann (b) 
boundary conditions and resonance modes in TM (c) and TE (d) cases, and these are 
specified by a mode index (±m, /) = (±8, 1). The eigenvalues and the resonance positions 
are k D{±8 ,i) = 12.2251 (a), k N ( ±8 ,i) = 9.6474 (b), k TM (±8,i) = 10.7845 - i0.02278 (c), and 
&te(±8,i) = 11.6295 — iO. 03993, where we take n = 2 in TM and TE cases. It is easily 
seen that the numbers of high intensity spots along the perimeter (angular direction) and 
along the radial direction are 2m and I, respectively. As expected from the above example, 
the eigenvalues and the resonance positions for one mode index (±m, /) have an order as 

kD(± m ,i) > Re[A; T E(± m ,i)] > Re[A; T M(±m,z)] > fcjv(±m,z)- 

We confirm numerically that all resonances in the TM case can be classified by the mode 
index (±m, /) just like the Dirichlet and Neumann cases. However, we find that in the TE 
case there are additional modes with a nonzero angular momentum m, absent in other cases, 
which can not be classified by the mode index (±m,Z). Figure 2 shows an example of the 
additional modes, (±m, I) = (±8, x) and kTE(±8,x) — 17.3507 — z2.4802. As expected from 
the figure and large absolute value of Im(A;), these modes are very leaky and originated from 
the existence of the Brewster angle in TE case on which rays can transmit without reflection. 
Therefore, these modes do not show any notable wave confinement by the dielectric interface. 

In order to investigate the distribution of resonance positions of the circular dielectric 
cavity, we obtain all solutions in the range of < Ke(kR) < 150 with R = 1 and n = 2 for 
both TM and TE cases. With the real part of the complex solutions obtained, we can check 
the level spacing distribution that is known to be Poisson and Wigner distribution for the 
integrable and chaotic billiards, respectively 6j. Although the circular dielectric cavity is an 
open system, it is still integrable, equivalent with the fact that all resonances can be specified 
by the mode index (m,l). We can, therefore, expect that the level spacing distribution is 
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Poissonic, and as shown in Fig. 3 the numerical calculation confirms this expectation. 

The properties of o pen ness in the circular dielectric cavity have been investigated by 
calculating the SPTD [13] which shows very different short time behaviors depending on 
the polarization. Especially, the exponential short time behavior appears in TE case and 
it has some relation to the existence of the Brewster angle. It is natural to relate this ray 
dynamical result with the distribution of imaginary values of the resonance positions since 
the imaginary value is relevant to the wave confinement by the dielectric cavity. 

The distributions of imaginary values of resonance positions for both TM and TE cases 
are shown in Fig. 4 (a) and (b), respectively. The substantial difference between both cases 
is that the distribution for the TM case is bounded, i.e., has a minimum imaginary value 



— 7Af(n)/2 while the result for the TE case is unbounded |18| . This result can be easily 



understood from details of reflection coefficients B,TM{n,6) and RTE{n,9), 9 is the incident 



angle, determined by the Fresnel equations 19]. Consider an initial ray with a fixed incident 



angle 9 which is invariant in the circular boundary. The survival probability of the ray would 

decay as I(t) = e" 1 ^ 1 ' 6 ' 1 where the decay rate is given by 

/ m In R TM (TE){n, 9) 
lM = 2^ ' (5) 

and the time t is scaled as the length of the ray trajectory. On the other hand, the intensity 

of a resonance decays as 

_ | e -*">*|2 = | e -i(Re(fc)+dm(fc))t|2 = e 2lm(fc)t ^ 

Therefore, the decay rate of a resonance mode (to, I) is given by 

l( m ,i) = -2Im(fc( m ,Q). (7) 

The ray dynamical (Eq.(jS])) and resonance mode (Eq.([7|)) decay rates have the same physical 

meaning due to the invariance of incident angle 9 in the circular boundary case. Therefore, 

if 7(n,6 l ) in Eq.flS]) is bounded, we can say that Im(A;( m n) is also bounded in TM case. We 

note that -Rtm(^, #) has its minimum at 9 = 0, i.e., the case of normal incident rays or 

bouncing ball trajectories. Then, the minimum Im (A>tm) is — 7A/(n)/2 and 

\nR TM (n,9 = 0) n + 1 
1m{u) = = In — . (8) 

The minimum imaginary value for the n = 2 case is — (ln3)/2 ~ —0.55 which is consistent 
with Fig. 4 (a). The above expression for the minimum imaginary value shows very good 



agreement with numerical results for various refractive indices, n = 2,4,6,8,10, as shown 
in Fig 5 (a). On the other hand, in TE case there is the Brewster angle 9b = arctan(l/n), 
and the ray incident with 9g transmits without reflection, i.e., Rte^i 9b) = 0. Then the 
decay rate 7(71, 9) in Eq.(jSD becomes infinity at 9 = 9b, and the imaginary value of resonance 
position in TE case can have very large absolute value of lm(k) like the additional mode 
(±m, x) mentioned before. This fact is consistent with the numerical result in Fig. 4 (b). 

We emphasize that the distributions of imaginary values of resonance positions, excepting 
the very high-Q resonances with only tunneling decay, can explain the behavior of the SPTD 
[13I ] . The resonance modes near the peak in Fig. 4 contribute to the short-time exponential 
decay in the SPTD, and the high-Q modes distributed just below the tunneling regime near 
zero in Fig. 4 do to the long-time algebraic behavior in the SPTD. In TM case, the fact 
that two mode-classes are smoothly connected, is consistent with the gradual transition 
from exponential to algebraic decays. On the other hand, in TE case, the distribution show 
an abrupt change at the peak, i.e., the two mode-classes are almost isolated, and this is 
responsible to the clear transition point from exponential to algebraic decays in SPTD. 



In order to understand who 
k( m ,i) with m = 40 in Fig. 5 (b) 



e shape of the distributions, we plot the resonance positions 



17 



181 ] . Black rectangle and red diamond represent resonance 
positions of TM and TE cases, respectively. It is known that the incident angle 9 of waves 
in a resonance mode (m, I) can be estimated by a semiclassical relation [?| 

sin0 = m/Re(A; (m j)). (9) 

From this relation we can obtain some Re(fc) values corresponding to the critical angle 9 C 
for total internal reflection and the Brewster angle 9b, which are indicated by black and red 
arrows, respectively, in Fig. 5 (b). If the incident angles of resonance modes calculated from 
Eq. ([9]) are greater than the critical angle 9 C , waves inside are very well confined by total 
internal reflection and then the resonance mode has very small |Im(A;)| as shown on the left 
of the black arrow. These high Q modes explain the peak near |Im(/c)| = in Fig. 4 (a) and 
(b), and the fraction of these high Q modes would be (1 — 1/n) in the semiclassical limit. As 
Re(/c) increases beyond the black arrow, lm(k) converges gradually to the minimum value 
— 7a//2 in the TM case, which corresponds to the bounded distribution near the minimum 
value in Fig. 4 (a). In the TE case, near the red arrow corresponding to 9b, Im(fc) shows 
very low value and then converges to —7^/2, which explains the unbounded distribution 
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and the peak at — 7m /2 in Fig. 4 (b). The same limit value — 7m /2 can be expected from 
the fact that Rtm^^O = 0) = RTE(n,9 = 0). The green and blue lines are obtained by 
calculating the relation between Re(fc) and lm(k) from Eq. fl^D, ([7]), and These explain 
well the resonance positions on the right range of the black arrow in Fig. 5 (b). 

Under the assumption that sin 9 values of resonance modes distribute uniformly such that 
P(sin#) = constant, we can obtain the distribution of Im(/c) of resonance modes as 

P(I,„ W ) oc («)-• (10) 
a{smv) 

Using Eq. (jSJ) and (J7]), we can calculate the distributions for both TM and TE cases, and 
the results are the red lines in Fig. 4 (a) and (b) which are in a good agreement with the 
histograms near the minimum imaginary value. 

Now we discuss the behavior of resonance positions at the small opening limit, i.e., 
n —>■ oo. Since the imaginary value of a resonance position kr± m n means the decay rate 
as shown in Eq. (JTj), we can safely take the limiting resonance position as a real value ki,. 
We first focus on the TM case. When n is very large, the solution of the boundary condition 
for TM polarization shown in Eq. (j4j) would be k = ki + 5 where 5 would be a small complex 
number and be zero at n — > oo. We can then expand the boundary equation Eq.(j4]) around 
ki- The Bessel functions can be written as 

J m {k L + 5) ~ J m {k L ) + J'Jk L ) ■ 5 + ■ ■ ■ , (11) 
= J m (k L ) + [-J m+1 (k L ) + —J m (k L )} ■ 8 H , 



where J'Jk L ) - 



<//,■ 



k=k L ) 



and the ratio of Hankel functions becomes 
H m (k L /n) 



(12) 



by the approximation by tangents of Bessel functions [15l . Il6j . The coefficients A an B is 
given by 

A - 2(m-l) ( m ] ' [ } 

B = h e lLfm-lr 1 )m-l/2 (M) 

2 V 2 ; V m(m- 1) ; V ; 

Then, the expanded equation of the boundary condition has small quantities, 1/n and S, 
and in the lowest order (1/n = and 6 = 0), the equation becomes 

J m -i(k L )=0. (15) 



This is nothing but the Dirichlet boundary condition (see Eq.(j4])), which means that at 
the small opening limit the resonance positions with order m of TM case approach the 
eigenvalues with order m — 1 of Dirichlet problem, i.e., 

lhn^ kTM(±m,l) — ^D(±(m-1),0- (16) 

From the equation for the next order, we know how the resonance positions converge to the 
limiting value k c . The result is 

5 TM = -An~ 2 -iBn- 2m . (17) 

Note that the imaginary part has an exponent depending on the angular quantum number 
m, while the exponent of the real part does not depend on the mode indices. 

For TE polarization, we can obtain similar results. In the lowest order, the TE boundary 
condition becomes 

J m {k L ) = 0. (18) 

This implies 

Note that at the small opening limit the resonance positions with order m of TE case 
approach the eigenvalues with the same order m of Dirichlet problem. From the equation 
for the next order, we get 

5 TE = -^ n ~ 2 - % {^L fBn^ m+1 \ (20) 
m m 

From the comparison between Im(5TAi) and Im(5^£;), it is clear that the TE boundary 
condition is more effective than the TM boundary condition in confining waves inside the 
circular cavity when the mode incident angle are greater than 6 C . 

In order to numerically confirm the above results, we trace two resonances of mode indices 
(m,l) = (8, 1) and (8,40) with increasing n for both TM and TE cases. As expected from 
the analytical results, we obtain 



lim k TM (8 ,i) 


= k D(7tl) ~ 11.0864, 


(21) 


lim £>TM(8,40) 


= k D(7A0) ~ 135.6942, 


(22) 


lim k T E(8,i) 


= k D{8 ,x) ~ 12.2251, 


(23) 


lim k T E(8,40) 


= A; D(8i4 o) - 137.2123. 


(24) 
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In Fig. 6 and Fig. 7, we show the n-dependence of 5tm and 5te- The results of Eq. flTT]) 
and fl20l) . denoted by solid lines, explains well the small opening limit behavior of resonance 
positions. It is noted that there is a transition point at n = n c ~ 17 in the traces of both 
&tm(8,40) an d /ctb(8,40) where the incident angle of the resonance modes(see Eq.Q) becomes 
the critical angle 8 C for total internal reflection. Therefore, the transition point indicates 
the change of leakage mechanism, from refractive to tunneling leakage of waves. For the 
refractive leakage range (n < n c ) , the trace of TM resonance modes shows 1/n dependence 

In summary, we have studied the distribution of resonance positions fc( m ,o hi the circular 
dielectric cavity, and the behavior of k/ m n at the small opening limit for both TM and 
TE polarizations. The distributions of lva.(kTu(m,i)) an d "k&(kTE(m,t)) are consistent with 
corresponding ray dynamical SPTD behaviors, and the small opening limits (n — > oo) 
of kj>M{m,i) an d hrE{m,i) approach different eigenvalues of Dirichlet problem, kD{ m -x,i) an d 
kD(m,i), respectively. The result implies that the small opening limit of some open system 
does not directly match to the corresponding closed system as shown in the TM case. 
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Figure Captions 



FIG. 1: The intensity patterns of a mode (m, I) = (8, 1) for various boundary conditions; 
(a) Dirichlet, (b) Neumann, (c) TM, and (d) TE boundary conditions. 
FIG. 2: The intensity pattern of the additional resonance mode (±m, I) = (±8, x) in 
the TE case. This corresponds to a very leaky mode showing wave escaping through the 
Brewster angle 9 b- 

FIG. 3: The level spacing distributions, (a) TM case, (b) TE case. These show good 
agreements with Poisson distribution (the solid line). 

FIG. 4: The imaginary value distributions of resonance positions when n = 2; (a) TM case 
and (b) TE case. Red lines are the results of a semiclassical analysis, Eq. CTUj) . 
FIG. 5: (a) The minimum imaginary values of the distributions for n = 2,4, 6, 8, 10 in TM 
case are denoted by solid dots. The solid line represents — 7m(b, 9 = 0)/2 and jm{ti, = 0) 
(see Eq. (JSj)) is the decay rate of bouncing ball trajectory, (b) The resonance positions with 
a fixed angular quantum number, m = 40, and n = 2 for both TM (rectangles) and TE 
(diamonds) cases. The left and right arrows represent Re(fc) values where the incident angle 
of resonance modes (see Eq. (jUj)) becomes to the critical angle 8 C and the Brewster angle 6 b, 
respectively. The resonance positions with incident angle greater than 8 C are well explained 
by results of the semiclassical analysis described by green (TM) and blue (TE) solid lines. 
Red cross represents an additional resonance mode which is originated from the existence 
of the Brewster angle. 

FIG. 6: The approaching behaviors of real part of the resonance positions corresponding 
to mode indices (m, I) = (8, 1) and (8,40) with increasing n (the small opening limit). The 
solid lines are the real parts (first term) of Eq. (ITT|) and (|20|) showing n~ 2 behavior. FIG. 
7: The approaching behaviors of imaginary part of the resonance positions corresponding 
to mode indices (m, I) = (8, 1) and (8,40) with increasing n (the small opening limit). The 
solid lines are the imaginary part (second term) of Eq.( fl~7l) and ( !20l showing n~ w (TM) and 
n- 18 (TE). 
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